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Abstract 

We develop a Melnikov method for volume-preserving maps with codimension 
one invariant manifolds. The Melnikov function is shown to be related to the flux 
of the perturbation through the unperturbed invariant surface. As an example, we 
compute the Melnikov function for a perturbation of a three-dimensional map that 
has a heteroclinic connection between a pair of invariant circles. The intersection 
curves of the manifolds are shown to undergo bifurcations in homology. 
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1 Introduction 

Volume-preserving maps on provide an interesting and nontrivial class of dynamical 
systems and give perhaps the simplest, natural generalization to higher dimensions of the 
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much-studied class of area-preserving maps. They also arise in a number of applications 
such as the study of the motion of Lagrangian tracers in incompressible fluids or of the 
structure of magnetic field lines [|l], Experimental methods have only recently been 
developed that allow the visualization of particle trajectories in three-dimensional fluids 
0, . The infinite dimensional group of volume-preserving diffeomorphisms is also at the 
core of the ambitious program to reformulate hydrodynamics 0. 

While some of the results for area-preserving maps generalize to the volume-preserving 
case, the study of transport in such systems is still in its infancy H, ||. The theory 
of transport is based on dividing phase space into regions separated by partial barriers 
through which flux can be measured. For the area-preserving case, the natural partial 
barriers are formed from the stable and unstable manifolds of periodic orbits or cantori 
@, |lOl- Primary intersections can be used to form resonance zones |1TT| , |12| — regions 



of phase space that are bounded by alternating stable and unstable segments joined at 
primary intersection points. Because the intersection points are primary, a resonance zone 
is bounded by a Jordan curve and has an exit and an entry set |TB[. The area of each of 
these sets is the geometric flux, the area leaving the resonance zone each iteration of the 
map. The images of the exit and entry sets and their intersections completely define the 



transport properties of the resonance zone [14 



Thus the beginning of a generalization of this theory to higher dimensions is the study 
the intersections of codimension-one stable and unstable manifolds for volume-preserving 
maps. 

As is well-known, a transversal intersection of stable and unstable manifolds is associ- 
ated with the onset of chaos, giving rise to the construction of Smale horseshoes. A widely 
used technique for detecting such intersections is the Melnikov method. Given a system 
with a pair of saddles, and a heteroclinic or saddle connection between them, the Melnikov 
function computes rate at which of change the distance between the manifolds changes 
with a perturbation. The integral of the Melnikov function between two neighboring 
primary intersection points is the first order term in the geometric flux p!6[| . 

Most applications of the method are for two-dimensional maps and flows |T^, 0, 
20|, though a Melnikov method for a three-dimensional incompressible flow was developed 
in In this latter case the perturbation may be periodically time dependent, and the 
Poincare map of the system is assumed to have a hyperbolic invariant curve, with two- 
dimensional manifolds. 

For the case of maps, the analogue of Melnikov integral is an infinite sum whose domain 
is the unperturbed connection. As usual, a simple zero of this function corresponds to a 
transverse intersection of the manifolds for the perturbed map. We developed a Melnikov 
method for three-dimensional maps in to study and classify intersections of stable 
and unstable manifolds for fixed points. 

In this paper we generalize this method to the problem of detecting heteroclinic orbits 
between a pair of normally hyperbolic invariant sets in volume-preserving maps on M". 
Our application is to the case of invariant circles for a three-dimensional map. 

To obtain a Melnikov function, we must define an appropriate measure of the distance 
between a manifold and its perturbation. In p2| we used the cross product of pair of 
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tangent vectors fields to obtain this distance. Different versions of Melnikov's method 
have used other ways of measuring the sphtting between the unperturbed separatrix 
and the perturbed one, though naturally only the normal distance is well-defined in the 
CO dimension-one case. This is appropriately measured using smadapted normal vector field 
or differential form. An adapted normal is a normal field to the saddle connection that is 
invariant under the dynamics. If the map is integrable, then the gradient of an integral 
can be used to construct the adapted normal, but the concept applies more generally to 
nonintegrable systems. 

We use the Melnikov function to construct a flux-form, an {n — l)-form whose inte- 
gral over a fundamental domain on the connection measures the first order flux through 
the connection. The fundamental domain is an annulus that generates the entire mani- 
fold upon iteration. Since the map is volume-preserving, the net (algebraic) flux always 
vanishes, but the one-way (geometric) flux gives a measure of the transport. 

In |^2| we introduced a family of volume-preserving maps that have a saddle connection 
between a pair of fixed points. This family is obtained from a family of planar twist maps 
with a saddle connection ||2^. This family can be modified so that it has a pair of 
invariant circles with a corresponding of saddle connection. We perturb this family by 
composing it with a near-identity, volume-preserving map, thus producing examples of 
volume-preserving maps with transverse heteroclinic orbits. 

We study the curves of zeros of the Melnikov function on a fundamental domain of the 
unperturbed manifold. Using the map to identify the two boundaries, the fundamental 
domain becomes a torus. Thus the zeros of the Melnikov function can be characterized 
by their homology on this torus. We show that as the parameters of the map are varied, 
the homology of these curves undergoes bifurcations, and that these bifurcations strongly 
influence the geometric flux. 



2 Basic definitions and properties 

Suppose /o : M" ^ M" is a diffeomorphism on n-dimensional Euclidean space. A smooth 
perturbation of /o is a family of functions fs = f{-,£) such that /(■, 0) = /o and /(x, e) is 
smooth in both variables. We now define a vector field on R" that will be used to measure 
the motion of an invariant manifold. 

Definition (Perturbation vector field). Given a perturbation f^ of /o, define the 
vector field for any point x by 

(1) 

Perturbation vector fields have some special properties. First, it is easy to see that 
is independent of /q. Second, if one regards as a time dependent vector field (where 
time is e), then y{e) = fe{x) is the solution of the initial value problem 

^ ^ , y(0) = /o(x) . 



d_ 

We 



fe{y) 
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Thus if we let = ft ° f^^, then F represents the flow of the nonautonomous vector 
field |2^, Thm. 2.2.23]. Since F is volume-preserving, the vector field has zero 
divergence with respect to Q [24, Thm. 2.2.24]. 

It is often convenient to define a perturbed family by composing /q with an e dependent 
perturbation: 

A = {td + ePe) O /o = /o + ePe o /o . (2) 

In this case, the vector field is the first order approximation to the perturbation, Xq{x) = 
Poix). 



2.1 Invariant manifolds 

Suppose the family has a family of invariant manifolds > M". In this paper, we 

will assume that We is a codimension-one surface. Our goal is to understand, at least to 
first order, the relation between the perturbation vector field and the way these invariant 
manifolds evolve with e. Later, we will restrict ourselves to the case in which We consists 
of pieces of stable and unstable manifolds of some invariant set. When We is a smooth 
graph over Wo, we can define a map that is adapted to the e parameterization: 



Definition (Adapted deformation). A map : Wo x (—eo,£o) 
We, if there is an Eq > such that 



i"-', is adapted to 



4>{-,e) is a diffeomorphism 0e '■ Wo 
0(-,O) = idyvo- 



We , V£ G i-eo,eo). 



There is quite a bit of freedom in the choice of (j); however, only the normal behavior 
is important for our application, since it measures the actual motion of We with e, and 
this is unique: 

Proposition 1. Suppose <^e andcl);; are two adapted deformations for a family of invariant 
manifolds We- Then 

-4>s{x) - T-Mx) e T^Wo. 



de 



de 



£=0 



Proof. Since both 0e(x) and (pe^x) are points in We, the curve Ce(x) = (t)~^{(t)e{x)) is a 
curve in Wq parameterized by e, and Cq{x) = x. Thus its derivative at zero is a tangent 
vector to Wo: 

|-Ce G T,.Wo. 

Thus 



§-/e{Cs{x)) 



e=0 



§-^Ux) 



+ DMx) 



e=0 



e=0 



^k{x) 



e=0 



Since D(f)o{x) = I, this gives the promised result. 



□ 
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We will use this proposition to compute the Melnikov function in For this we 
need to measure rate of change of an invariant manifold with respect to a perturbation — 
changes in the tangent direction are unimportant. In order to measure the change in the 
normal direction, we introduce the concepts of adapted normal vector fields and adapted 
forms. 

2.2 Adapted normals 

We will measure the splitting by using a normal to the invariant manifold Wo, that 
throughout this paper will be a codimension-one submanifold, i.e. a surface. To be 
useful, the normal field should evolve in a precise way under the unperturbed map or, as 
we say, be "adapted" to the dynamics. 

First we recall some notation. Let f be a vector field and / a diffeomorphism. The 
pull-back of V under / is = {Df{x))~^v{f{x)). Similarly, the pull-back of a 

/c-form a; is {f*uj)x{vi,V2, ■■■Vk) = ^f(x){Df{x)vi, . . . ,Df{x)vk). Finally the inner product 
of a vector field Y with a fc-form uj is defined as the {k — l)-form iyu: = uj{Y, -,...,■). 

Definition (Adapted normal field). Suppose that / : M" — M" is a diffeomorphism 
with an invariant surface W, and there is given an inner product (, ) for vectors on M". 
An adapted normal field is a smooth function : W — » M" such that 

• r]{x) 7^ for all x G W. 

• ri{x) is normal to the surface for all x G W, that is t]{x) G T^.W"'". 

• For all vector fields y : W — * M" we have that 

f*{r^,Y) = {r^J*Y). (3) 

The geometry is shown in Fig. |l]. Note that the puUback of a scalar function (7 : ^ M 
is f*g{x) = g{f{x)); thus if we define Z{x) = f*Y{x), @ is equivalent to 

{v{x),Zix)) = {r]{f{x)),Dfix)Zix)) 

Adapted normals can be thought of as a generalization of the gradient normal that 
one gets from a first integral. Recall that the gradient of a smooth function J : M" M 
is the unique vector field, VJ, such that for all vector fields Y on R", 

iydJ = dJ{Y) = {VJ,Y). (4) 

If / has a nondegenerate first integral, J = Jo/, then (^) implies that f*{'VJ,Y) = 
(V(J o /), f*Y). Therefore, if the diffeomorphism / has a first integral J, then VJ is an 
adapted vector field, provided it doesn't vanish on W. 

If we are using the standard inner product on W\ then we can characterize adapted 
normals more concretely. 
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fYix) 




'X 



fix) 



Dfix) 



,^Y(f(x)) 



Figure 1: is an adapted normal if the two rectangles shown have the same area. 



Proposition 2. Let 77 : W ^ R" be a smooth function defined on the invariant surface 
W, and suppose {u,v) = ■ v is the standard inner product on M". Then t] satisfies (Qj 
for all vector fields F : W — > M" if and only if, for all x 

Df{xMf{x)) = r^{x). 

In the general case W is not defined as the level set of an invariant, and it is not easy 
to show that an adapted normal field exists. However, when the map is volume-preserving 
and we are given an appropriate parameterization of the invariant surface, an adapted 
normal vector field can easily be constructed. 

Lemma 3. Suppose that / : — > M" preserves the volume-form Vt, and has a smooth 
invariant surface W. Suppose k : M"^^ — > M" a nondegenerate parameterization ofW 
with the property 

f{k{u)) = k{u + 5) , (5) 
for a constant 5 G M"^"'^. Then the vector field 77 : W ^ M" restricted to W defined by 

{r], ■) = du^k, du^k, . . . , 9„„_,/c) (6) 

is an adapted vector field on W. 

Proof. The nondegeneracy of the parameterization implies that 77 7^ on the surface W. 
Condition (^) implies that 



D^k{u + 6)=D^{fo k) {u) = DJ {k{u)) D^k{u) ; 
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therefore, since = we have 



r{r^{k),Y{k)) = M{k)),Y{f{k))) 

= {Y{f{k)), duMu + 5), . . . , du„_Mu + 5)) 
= {Y{f{k)), Df{k)d^^k{u), Df{k)du^^Mu)) 
= {Df{x)-'Y{f{k)), duMn), . . . , du„_,k{u)) 
= {r^{k),rY)) . □ 



2.3 Adapted one- forms 

An alternative concept to that of adapted vectors are adapted forms. The advantage of 
one versus the other approach is mainly a question of taste, though differential forms can 
be used without assuming an inner product. We now define an "adapted one-form." 

Definition (Adapted one- form). Suppose that / : M" — > zs a diffeomorphism. An 
adapted one-form on an invariant surface W is a smooth function v : Ty^W — > M such 



• Vx is nondegenerate for all x G W. 

• u^{v) = for all v G T,W. 

• f*v = v 

Note that when v is an adapted one- form then for each x G W, ker(z/a,) = Tj.W, 
but that since v is nondegenerate, it will not be zero for vectors that are not tangent to 
W. As before, we note that if W is given as the level surface of an invariant function 
J, i.e. if J{f{x)) = J{x), then an adapted one- form is easy to obtain: the one- form 
dJ is adapted provided only that J has no critical points on W. This follows because 



Given an inner product, (■, ■) we can always associate a unique vector field, t] with a 
form u, through ixi^ = t^iX) = {ri,X). Here : W — is a smooth function. It is 
easy to see that, if the rj is adapted, then z/ is also adapted. Conversely, given an adapted 
one-form, we can find an adapted normal field, through the same relation. 

Proposition 4. Let rj and v he related through 



Corollary 5. Let / : ^ preserve a volume-form VL, and W be a smooth invariant 
surface. Suppose k : M"~^ ^ is a nondegenerate parameterization of W such that 
f{k{u, v)) = k[u + 5) for constant 5 G M"^"*^. Then 



that 



f*dJ = d{Jof) = dJ. 



Then rj is an adapted normal if and only if v is an adapted form. 
Using this with Lem. |^ implies immediately. 



(7) 



^(■) duik, . . . , du„_ik) 




is an adapted one-form on W. 
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2.4 Example 

Let / : M'^ — s> be the two parameter family of diffeomorphisms 

(e'^ {x cos 9 — y sin 9) \ 
e'^ {x sin 9 + y cos 9) , 

where r and 9 are constants. It is easy to see that / preserves the standard volume-form 
Q = dx A dy A dz. In addition, / has an invariant surface given by 

W = { (x, y, z) e \ {0} : + = 4z sinh(2r) } ; 

however, the "obvious" function + — 42;sinh(2r) is not invariant. Instead, we pa- 
rameterize W with the function /c : ^ given by 

(e"'^cosf \ 
e" ^ sin V 
(4sinh2r)-ie2"^ / 

The function k is nondegenerate and satisfies 

f{k{u,v)) = k{u + l,v + 9). 

so that Lem. ^ and Lem. ^ apply. Using (^, we find that u = —2Tz{xdx + ydy — 
2sinh.{2T)dz) an adapted form on W. In other words /*z/ = u and kei Up = TpW, as can 
be explicitly verified. Also using (^, it is possible to show that 

/ -2txz \ 
r]{x,y,z) = -2Tyz 

y Atz sinh 2r J 

is an adapted normal field on W. In other words, it satisfies Df{xyr]{f{x)) = ri{x), for 
each X G W. Note that u and rj are related through ixi^ = {r],X). 



3 Melnikov Function 

Suppose that the diffeomorphism /q has two normally hyperbolic invariant sets p and q, 
and a codimension-one surface W = = W^{q) that is a saddle connection between 

them. Upon perturbation, suppose that the corresponding invariant sets Pe and of fs 
have a stable manifold W| and unstable manifold W". Then in our notation, the classical 
Melnikov function is the smooth function M^, : W ^ R on the saddle connection W 
defined by 

d 



de 



e=0 



(9) 
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for a given adapted form u on W, and a given pair of adapted perturbations 0| and 0^ 
corresponding to the stable and unstable manifolds, respectively. Thus measures the 
relative "velocity" of the manifolds as a function of e. While M^, appears to depend on 
the choice of adapted perturbations, we will show that it does not. 

A function similar to M^, was used in ||2^, to study the topology of heteroclinic con- 
nections of fixed points. Our purpose is to apply the method to the case of invariant 
circles, as illustrated in Fig. ||. 




Figure 2: Two normally hyperbolic invariant circles Ci and C2 with a saddle connection W. 



3.1 The fundamental iterative relation 

The fundamental relation used in deriving the Melnikov function is an iteration formula 
obtained by combining the definition of adapted one-form and Prop. |1]. 

Theorem 6. Suppose is a family of diffeomorphisms with invariant surfaces Ws- Let 
V he an adapted one-form and 0^ he an adapted deformation on Wo- Define /i : W ^ M 
hy 

Then 

/^-/io/-i = zy(Xo), (10) 

where Xq is the perturhation vectorfield, Moreover, if (pe is another adapted deforma- 
tion and jl is defined similarly to fi, then fJ^ = fi- 
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Proof. By Prop. the difference between the derivatives of two adapted diffeomorphisms 
is tangent to W, and by Defn. |2]^, z/ vanishes on any tangent vector. Thus fj, = fl. To 
compute the second relation, use /qZ/ = z/ to find 

= foiyxide(l)eix)) 

= Z^/(,(x) {Dfo{x)de(l)e{x)) 

where we suppress the e = expressions for simphcity. Using (|I|) Xq = defe{fQ^{x)), we 
have /i o f-^ = [dsfs{(f>s{fo^{x)))) - z/^(Xo), and therefore 

Noting that 0^ = fe{(f>e{fQ^{x)) is also an adapted diffeomorphism, we see that the first 
term vanishes by Prop. |l|. □ 
Equation (|10D gives us a recursive formula to compute the normal component of the 
change in the manifold W^. 

Corollary 7. Under the assumptions of Thm. for all n 

n-l 

li = liof^- + Y,^{Xo)of-K 

k=0 

In addition, if lim fi o /o~"(x) = 0, then 

n— >oo 

oo oo 

/x(x) = K^o) o f,~' = J2 ^ iiu'r^o) . (11) 

fc=0 fc=0 

These statements can be directly transcribed for adapted normals using Prop. ||. 
3.2 Transversal intersections 

According to Thm. ^ and Cor. we can compute the Melnikov function (H) in terms of 
the first order perturbation vector field Xq. 

Proposition 8. Suppose f has a codimension-one saddle connection W between two 
normally hyperbolic invariant sets p and q. Assume that for all x G pUq, the perturbation 
vector field Xq{x) = 0. Let v be an adapted form and rj the corresponding adapted normal 
defined on W. Define the Melnikov function by 

oo oo 

M,= J2 K^o) o /o^ = Yl (V, Xo) o f^ . (12) 

= — oo fc = — OO 

Then if a point xo & W is a nondegenerate zero of Mi,, the stable and unstable manifolds 
IV"(g, /e) and W^{p, fs) intersect transversally near xq for e small enough. 
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Proof. For each point x in the saddle connection W, there is a neighborhood Aq C W, such 
that all the iterates /'^(Aq) are disjoint. Moreover, since p and q are normally hyperbolic, 
the stable manifold theorem implies that there is an eo > such that there exist adapted 
deformations 0" : ATq x (-^o, ^o) fe), and 0^ : A/'o x {-sq, Eq) W'{p, /J. 

Consider first the unstable part. Let V = [JkLo fo^i-^o)- Clearly V is a immersed 
manifold. Moreover, we can extend the domain of to all of V, by defining 

provided that x G /q"^(A/o). It is clear that for each e G (— ^o^^^o) and x G V, we have 
that G WiqJ,). 

For each x, we are interested in estimating (j){x, e) to first order in e. Using 0" in 
Cor. 1^ gives ( |Tl| ) providing fi'^ o /""(a;) 0. This is the case because 
so that de4>^ is bounded, and J^f"(^x) ~^ since it is an adapted form. 

Similar analysis applies to the stable adapted deformation, and again Cor. |^ applies, 
though we iterate in the opposite direction, to obtain = — Xlfcli '^(-^o) ° /o • According 
to (^, the difference between /x" and /i^ gives the Melnikov function, which yields (|T2p. 

Following a standard Melnikov argument based on the implicit function theorem , 
we conclude that if xq is a nondegenerate zero My then near xq, the two manifolds W^"(/e) 
and W^{fe) intersect transversely. □ 



4 Flux 

The flux across a surface is the volume that crosses the surface each iterate of a map; it is 
an important measure of transport. Recall that for area-preserving maps, the Melnikov 
function is a measure of the distance between the stable and unstable manifolds, and that 
its integral between two successive zeros is the geometric flux that crosses the "separatrix" 



each iteration of the map |jT5|, [l^]. The outgoing flux is exactly balanced by an ingoing 
flux, so that the net, or algebraic, flux crossing the separatrix is zero. 

Here we will obtain an analogous formula for volume-preserving maps (see also [|]). 
We start by constructing a flux form on an invariant set. We will see that the algebraic 
flux crossing the separatrix is zero. This implies, for example that the Melnikov function 
has zeros in the separatrix. 



4.1 Flux Form 

It is well known that a volume-preserving map with an invariant J can be restricted to 
a measure preserving map on any surface J = c on which VJ is nonzero. That is, the 
form uj = \VJ\^'^i^jQ is an invariant n — 1 form for the map /| We show here that 
a similar preserved measure also exists if we can flnd an adapted normal for an invariant 
surface W. We will then use this to construct a flux form on W. 

As usual, we assume that / is a diffeomorphism with an invariant volume-form Q, W 
is an invariant codimension-one hypersurface, and ( , ) is an inner product on M". 
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Proposition 9. Suppose rj is an adapted normal field on W. Then 

{ri, rj) 

is a nondegenerate {n — l)-form on W that is invariant under the restricted map /|yy : 
W ^ W. 

Proof. It is clear that cj^ is a nondegenerate [n — l)-form on W. We need to show that 
/qC(J^ = for vectors in TW. With some manipulations we have 

f*UJjj — UJrj 

J \'h III \'h III 

iyil. 

where we define the vector field: 



if*^f*n 




f*{v,v) 


{v,v) 


i j.^n 




[r], f*r]) 






V 



Since {t], w) = on W for each point x G W, and rj defines the normal direction, then 
v{x) e TxW. Since v is tangent to W, the form iyfl restricted to W has to vanish. Thus, 
we conclude that {f*ujr^ — uJn)y^ = 0. □ 
From now on assume that W has an adapted vector field rj, and the Melnikov sum 
(1^) exists. We then define 

Definition (Flux Form). The flux form $ = Mr^uJ^^ is an {n — l)-form on W. 

Note that $ might be degenerate (it has zeros), and since the space of (n — l)-forms on 
the {n — l)-dimensional manifold W is one-dimensional, it might not be unique. However, 
this is not the case: 

Lemma 10. // the Melnikov function exists, then the form $ is independent of the choice 

of 7]. 

Proof. The projection of Xq onto TW is the vector v = Xq — Note that 

ixo^ - {r],Xo)uJr, = iyVl. 

Since v G TW, we conclude that i^Q = 0, as a — l)-form in the surface W, and so we 
have 

{r],Xo)uJr, = ixo^ 

for all vectors in TW. This implies that the summand of MrjUJr], recall (|12|) , can be 
rewritten 

if^r{v,x,)uJ, = if^r^x,n, 

which is independent of rj. □ 



4 FLUX 



13 



Thus we have 

oo 

$ = M,u, = U'oT^Xo^ ■ (13) 

fe=— oo 

Since is nondegenerate, the degenerate points of the flux form correspond to zeros of 



the Melnikov function. As we will see in §4.3| , the integral of the flux form over a piece of 
W gives flux through that surface to first order in e. 
The form ix^^ has some interesting properties: 

Proposition 11. Let be a perturbation vector field Then the form ix^^ is exact. 



Proof. As we already noted, Thm 2.2.24 in implies that the divergence of X^ vanishes 



which by definition {dwQX^)^l = means that the Lie derivative vanishes as well. 

Since dQ = and Lx^^ = d{ix^^) +ixsdQ, this implies that d{ix^^) = 0. Thus the form 
is closed. Since ix^^ is globally defined in M", the form is exact. □ 
Using this result we can obtain an {n — 2)-form on W such that 

d(3 = ix,^ ■ (14) 

In this case we will say that (3 represents the perturbation on W. Using Lem. lU and (|T^, 
it is easy to see that if exists then the following [n — 2)-form is well-defined on W. 

oo 
fc=— oo 

Notice that a is invariant under / and is independent of rj. Using this we can see that: 
Proposition 12. The flux form $ is 

^ = da . 

Proof. This is a straightforward calculation using Lem. ITO and Prop. 0. □ 



4.2 Fundamental Domains 

Our goal in this section is to find a compact subset of the manifold — a fundamental 
domain — that generates the entire manifold under iteration by /q. We will integrate 
the fiux form over the fundamental domain to show that the algebraic fiux crossing the 
separatrix is zero. From this point on, we will concentrate on the case n = 3. To define 
the fundamental domain we start with the concept of a proper loop: 

Definition (Proper loop). Let f^-.M?^ be a diffeomorphism, and W a forward 
invariant surface. We say that a smooth Jordan loop 'j G W is a proper boundary in W 
if there 7 bounds a surface W-y C int(yV) which is a trapping region: 

f{d{Wj)) C mt(W^) . 
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Similarly a loop 7 is a proper loop for a backward invariant surface if it is a proper 
loop for the map /^^. 

It is important to notice that not all invariant surfaces admit proper boundaries. A 
trivial observation is: 

Proposition 13. If is a proper boundary in W, then f{'j) is also a proper boundary. 
In addition, W/(^) = /(W^). 

The situation that we have in mind relates to the structure of stable and unstable 
manifolds. Let a, b be compact, normally hyperbolic invariant sets of /, and W = W^{a) = 
W^{b) a saddle connection between them. A proper loop 7 C W is a submanifold of W 
that bounds a local submanifold that is a isolating neighborhood of a in W^{a). In other 
words 7 is proper if it bounds an open local submanifold, Wf^^i^a) = W-y, that maps inside 
itself. 

If 7 is proper, we can define the stable manifold starting at'y, denoted by = W^{a), 
as the closure in W^{a) of the local stable manifold bounded by 7. In the same way, for 
b, if we have a proper loop a for f~^, we define the unstable manifold up to a, denoted 
W^{b), as the interior of the local unstable manifold bounded by a. We will see below 
why it is convenient to use this slightly asymmetric definition. 

Given a proper loop we can define 

Definition (Fundamental domain). Let W be a forward invariant surface. An sub- 
manifold with boundary, V, is a fundamental domain of W if there exists some proper 
loop 7 in W, such that 

The fundamental domain is a manifold with the boundary 

9P = 7U/(7), 

see Fig. |^. An immediate consequence of the definition is that all the forward iterations 
of a fundamental domain are also fundamental and Vj^^) = f(V.y). It is easy to see that, 
if proper boundaries exist, then the forward invariant manifold can be decomposed as the 
disjoint union of fundamental domains. 

w = (>v\>v-,)u|J/^(p) . 

fc>0 

If the surface W is both forward and backward invariant, then this decomposition works 
in both directions. In such case we have 

^=(jf'iV) . 
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Figure 3: Fundamental domain V on a stable manifold of an invariant circle a bounded by a loop 7 and 
its image /(7). The second part of the figure shows the annular fundamental domain itself, together with 
an assigned orientation. Finally, if we identify the points on 7 with their images, then the fundamental 
domain is equivalent to a torus. 

4.3 Algebraic Flux 

Given a vector field X, the differential form ix^ represents the flux associated with 
X] that is given set of vectors f 1, f2, ...fn-i, ^{X,vi,V2, ■■■Vn-i) is the volume of the 
parallelepiped formed from these vectors, and thus measures the rate at which volume is 
swept out by X through the parallelepiped defined by vi, V2, 

According to (|13|), the form $ is the sum of ixo^ along an orbit on W. Thus $ 
evaluated at a point on a fundamental domain V measures the total flux of Xq along the 
orbit of that point. 

The algebraic flux through a surface is the integral of the flux over the surface. Since 
$ measures the flux along an orbit on W, the integral of $ over a fundamental domain is 
the algebraic flux through the entire surface W. 

Proposition 14. The algebraic flux through W is zero: jy^ix^P- = Jp^ = ^■ 

Proof. The fundamental domain V = is a. submanifold with boundary, such that 
dV-y = 7 U /(7), where 7 is closed curve that does not intersect /(7). If we give an 
orientation [V] to V, the induced orientation on the boundary satisfies [7] = — [/(7)], 
recall Fig. |^. Since ^ = da hj Prop. |T2], and a is invariant under /, Stokes's theorem, 
implies 
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A simple corollary is 

Corollary 15. The Melnikov function must have zeros on V. 

The importance of fundamental domains is that much of the information about the 
entire manifold can be found by looking only at these submanifolds. In particular, if we 
have a transversal intersection of two invariant surfaces, we can look at a pair of funda- 
mental domains and study primary intersections. In Fig. |^ we show a pair of fundamental 
domains of two stable and unstable manifolds that intersect transversally. 



Figure 4: A pair of fundamental domains for two different circles Ci and C2 that intersect transversally, 
forming three-dimensional lobes. 

The curves of zeros of the Melnikov function can be classified by their homology on V. 
To do this, we identify two boundaries of the fundamental annulus by identifying 7 with 
/(7). With this identification the fundamental annulus becomes a torus, as sketched in 
Fig. ^. Since the homology group of the torus is Z^, we can label the curves by a pair of 
integers (m, n) which represent the number of times the curves wrap around each circuit 
of the torus. For example when the identification is performed on Fig. ^ there are a pair 
of zero crossing curves with homology type (3, 1) — they move once around the annulus in 
three vertical transits. 




a 



b 



5 Examples 



In this section we construct a family of volume-preserving maps that have a saddle con- 
nection between a pair of invariant circles. We obtain this family by starting with an 
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area-preserving twist map that preserves an axis and extending it to a three-dimensional, 
volume-preserving map by composing it with a sheared rotation about that axis. The 
twist map is defined in such a way that it has a saddle connection between two fixed 
points, and so the resulting three-dimensional map has a pair of invariant circles with a 



two-dimensional connection. Examples similar to these were found by Lomeli and are 
closely related to those in p^ . 

We begin with an area-preserving map on in coordinates {z, r) that preserves the 
axis r = 0, and has a fixed point at some nonzero r = r*. For example, set 

(r', z') = G{r, z) = {h^^{r + h{z)) — z, h{z) + r — r*) , 

where r* G Z. Here we assume that : M ^ R is an increasing circle diffeomorphism of 
period 1, i.e., h{z + 1) = h{z) + 1. Moreover we can verify that det{DG) = 1, so that G 
is area-preserving. Finally 

G-\r, z) = {z- h{h-\z) - r), h~\z) + r* - r) , 

so that G is a diffeomorphism. 

It is easy to see that G{0, z) = (0, h{z) —r*) so that the 2;-axis is preserved. The map 
has fixed points at solutions of z = ^{h^^{z) + h{z)), with r = r* + z — h{z). In particular, 
any hyperbolic fixed point of h, z* = h{z*), yields a saddle fixed point (r*, z*) of G whose 
multipliers are A = h'{z*) and 1/A. Between every pair of such fixed points of h there is 
at least one other fixed point of G; it is typically elliptic. 

The map G is not necessarily integrable (in §^]T] we will choose an h that leads to an 
integrable map). However, G always has a pair of invariant curves: 

Wo = {iz,r) ■.r = r*} , 

Wi = {{z,r):r = h^\z)~h{z)+r*} . (15) 

These curves intersect at any fixed point z* of h. Thus they provide a saddle connection 
between points {zl,r*) and {z2,r*), where z* are neighboring fixed points of h. 

We show an example of the dynamics of G in Fig. |^ for the case that h{z) = z — 
^cos(27rz), and r* = 1.0. Here one can see the saddle connection at r = r* as well as 
chaotic dynamics in other regions of phase space. 

We can extend G to by introducing the cylindrical angle 6 and using the volume- 
form fl = dr A d9 A dz. Defining r to be the "symplectic" polar radius, 

r = ^ix^ + y^), (16) 

the cylindrical coordinates are 

(x, y, z) = P{r, 6, z) = ( v/2r cos 6, V2r sin 6, z) 
so that Q = dx Ady A dz. In terms of these coordinates the map becomes g = PoGo P^^: 



g{x, y, z) = {p{r, z)x, p{r, z)y, r + h{z) - r*) , 
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where p = ^Jr' jr is explicitly 



h i(r + h{z)) - z 

1 T ^ \} 



It was shown in |22| that if h is C"', then p is C^'^ so that, in this case, (7 is a 
diffeomorphism. 

The map becomes fully three-dimensional if we introduce dynamics in the angle Q. To 
do this, we compose the map with a rotation. Denote the rotation about the z-axis by 
angle by 

(cos — sin \ 
sin-?/' cos?/' . (18) 
01; 

Introducing a rotation angle T{r,z) that depends smoothly on {r,z), we define a diffeo- 
morphism / by 

f = 90Rr, (19) 

Note that since Rr and (7 both preserve fl, so does /. 
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The map / still has a rotational symmetry 

foR^ = R^of, (20) 
for any constant G M. This implies that when G has a saddle connection, so does /: 



Proposition 16. The surfaces ( 11^) are invariant under (\T^). In addition, Wo and Wi 
intersect on the invariant circles 

C{z*) = {{x,y,z) : z = z*,r = r*} , 

where z* is any fixed point of h. 

Every point on the circles C{z*) is fixed under g. The derivative of g at such points is 



Dg{x,y,z) 







v 



X 



y 



(21) 



where A = h'{z*). More generally, we can compute the derivative of /" on the invariant 
circles for the special case that the rotation angle is constant: 

Proposition 17. Suppose f is given by (|7^, that t is constant, z* = h{z*),and C{z*) is 
the corresponding invariant circle. Then for all {x,y,z) G C{z*) 



Dr{x,y,z) = R„ 



2,.v- - . y') 2ir(A-"-l)xi/ \ 
^{X-^^-l)xy ^(a;2 + A-V) 
(A2"-l)x (A2"-l)y 



V A"-i(A2-l) A"-i(A2-l) 



A" 



(22) 



where A = h'{z*). Moreover, if X > 1 {< 1) the invariant circle has a stable (unstable) 
manifold contained in Wi, and unstable (stable) manifold contained in Wo- 



Proof. Given the symmetry (POI), it is enough to check (E^) for points of the form 



^,0,z*). 



Since (^2]) reduces to ( ^1]) when n = 1, it is enough to verify the induction step 
Dg{V2^, 0, z*)Dg'\V^, 0, z*) = Dg''+\V2^, 0, z*) . 

The vector (0,0,1)* is an eigenvector of Dg{x,y, z*) with eigenvalue A. Since this is 
tangent to Wq, this shows that it is the stable manifold when A < 1. Similarly the vector 
{x,y,2r*jA^y is an eigenvector with eigenvalue A~^ that is tangent to Wi. The final 
eigenvector of Dg is {y, —x, 0)* which is tangent to C and has eigenvalue 1. □ 
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5.1 Integrable case 



In general the maps g and / are not integrable, even though they have a saddle connection. 
However, for a special choice of h there is an integral. This example is related to the 
work of Suris |p6i p7| on area-preserving integrable maps, but is distinct from the three- 



dimensional maps found in p8[ that have an invariant but which do not have a rotational 
symmetry. 

Let m(w) = 3^ be the Mobius transformation on R U {oo} with ad — be = 1. A 



cw+d 



circle map conjugate to m is obtained by defining w = tanvr^;, giving h 



hm{z) = — arctan (m(tan7r2;)) 



7r 



(23) 



This map can be written more explicitly as a circle map using trigonometric identities: 



hm{z) = z ^ — arctan 



h — c + {h + c) cos(27rz) + (a — d) sin(27rz) 
a + d — [a — d) cos(27rz) + (6 + c) sin(27r2;) 



(24) 




Figure 6: The circle map h^iz) and its inverse for (a, b, c, d) — (1, i, 2, 2). 

Requiring ad — bc = l, some useful properties of this family of circle maps follow easily 
from its conjugacy to the Mobius transformation 

• hl^ = h^t, for all t G Z; 



If |tr(m)| > 2 then hm has two fixed points G [0, 1) 
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• The fixed points have multipliers Dh{z^) = j^^^^^^r^^ > 0, where a = a + d. 

• Thus is unstable and z"^ is stable. 

For any hm, the resulting map G is integrable. To see this, we show how this map is 
related to the Suris example. First note that G can be rewritten as a second difference 
equation 

zt+i - 2zt + zt-i = h{zt) + h'^{zt) - 2zt = F{zt) , 

where = Zt+i — h{zt) + r* . Suris showed that this family is integrable when F is given 
by 

, , 1 r Asin(27rx) + B cos(27rx) + C sin(47rx) + D cos(4'7rx) 

F(z) = - arctan \ ^ ^ / ^ ^ r ^ / 

vr [1 - E - A cos(27rx) + B sin(27rx) - C cos(47rx) + D sin(47ra;) 

for any values of the parameters A, B, C, D, E. After some algebra one can see that our 
map has this form with A = b"^ — cP, B = (a — d){c — b), C = |(6 + c)^ — |(a — d^, 
D = {d-a){b + c), E= i(a2 + b^ + + d^). 
The map G with this h has the integral 

J{z, z) =(1 - E) cos(27r(2 - ^)) - A(cos(27rz) + cos(27r2')) + 

S(sin(27rz) + sin(27r/)) - C cos(27r(z + z)) + D sin(27r(^ + 2')) • 

For the examples, we will use the map (|1^) with h.^ given by (^) with 

m{w) = ' 25 

This corresponds to setting m to the hyperbolic rotation matrix a = d = cosh(ln(z/)), 
b = c = sinh(ln(i/)). In this case m* is given by replacing u with z/*; thus iteration of h is 
extremely easy. This was also the example used in [^. 

This gives a family of three-dimensional maps, /, with parameters i' and r. Setting 
r* = 1, there are invariant circles at {r*,z*) = (1, ±i). For this case the invariant has the 
form 

J(x, y, z) = 2u cos(27rr) + (1 — z/^) cos(2'7r2;) sin(27rr) . (26) 

The level sets corresponding to J = 2v give the invariant manifolds Wq, Wi, and the 
circles C{z*). The level sets of J are shown in Fig. |^ for the case v = 0.3. Though the 
level sets of J make it appear that (|, ±^) are also invariant; G(|, \) = (|, —j), so these 
points move downward. Moreover the curve {r = i} has image {r = h~^{z) — hy{z) + |}, 
and this latter curve has image again of r 



1 

2" 



5.2 Perturbed Map 

We break the invariant surfaces by choosing a perturbation of the form (Q). For the first 
example, we choose a composition of two simple perturbations: 



Pi(a;,i/,z) = ((l + i/2)(^*2_^2)^Q^Q) 
P2{x,y,z) = {Q,x{z*^-z^),Q) . 



2 2, .^ (27) 
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Each of these maps has a nilpotent Jacobian, which imphes that the maps id + ePi are 
volume-preserving for all e. The complete perturbation is then defined as 

id + ePe = {id + ePs) ° {id + ePi) . 
Substituting the perturbation into the computation (0) for the vector field gives 

Xo = P2 + Pi . 

Fundamental domains Vi on Wj are given by the annuli bounded by the circles 7^ = 
= 0} n Wi and their images, /(70) = {z = h{Q)} n Wo or /(71) = {z = /i-^(O)} n Wi, 
respectively. These can be projected onto (2;, 0) coordinates for visualization. Calculation 
of the Melnikov sum (JT^) is straightforward using the adapted form dJ associated with 
the invariant (PB]). 

We show several representative contour plots of Mdj in Fig. ^ In the figure, positive 
values of Mdj are shown as dashed lines and negative as dotted lines, while the zero 
contour is the solid line. For example, in the bottom-left panel {v = 0.275 and r = 0.325) 
there are two zero contours, corresponding to the unstable and stable manifolds crossing 
with opposite signatures (since the algebraic flux through the fundamental annulus is 
zero, the zero contours must come in pairs). Since the unperturbed map takes the circle 
71 to the circle /(71) shifting each point by r, the lower boundary of V can be identified 
with the upper boundary after shifting the latter to the right by r (we show this shift by 
the arrows in Fig. After this identification the fundamental annulus becomes a torus, 
and the zero contours correspond to a pair of circles that wrap once vertically. Thus these 
contours have homology type (1,0). The bottom right panel also has this homology type, 
though the curves are very close to a bifurcation. 
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v=0.275, t=0.325 



v=0.25, T=0.425 



Figure 8: Contours of the Melnikov function for Wo with the perturbation (27). Shown are four values 
of the parameters {f, r). Bounds for the figures are z — [0, h~^{Q)] and 9 = [0, 2tt]. The arrow at the top 
of each panel shows the translation by r. 



There are several such bifurcations in homology type of the zero contours as we vary 
the parameters. For example in the upper left panel, the homology type is (3, 1) — as each 
zero contour moves from the bottom to the top of V, it lags the maps translation of 6 by 
a full circuit in three vertical transits. In the top-right panel there are two zero contours 
with the homology (2, 1). To elucidate these changes in homology, we show a bifurcation 
diagram in the space of the parameters in Fig. ^. We have only found the three homology 
classes already mentioned. 

Also shown in Fig. ^ are contours of the geometric flux 

Flux = - / |<l>| 

as a function of u and r. The flux is largest when u and r are both small, and it appears 
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to get extremely small as u approaches one. Note that there is a "valley" in the flux 
contours near both homology bifurcations. 
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Figure 9: Contours of the geometric flux through Wq as a function of and t for the perturbation 
(p7|). The nine contours are at equally spaced levels ranging from a flux of 0.09 at the lower left to 0.01 
at the top. Also shown are bifurcation curves corresponding to the change in homology types of the zero 
contours of Mdj- 



Finally, we have also studied the perturbation 

P,{x,y,z) = {{1 + y'){z*'' - z^), 0,0) , 

P2ix,y,z) = {0,x'iz*'-z'),0) , (28) 
Ps{x,y,z) = iO,0,r-r*) . 

giving a perturbation vector field Xq = Pi + P2 + P3. We show the bifurcation diagram 
for the zero contours of M^j for Wq in Fig. 10. For this case there appear to be only two 
homology types, (1, 0) and (3, 1). Again there is a "valley" in the flux near the bifurcation 
curve. 

We have also computed the Melnikov function for the second invariant set, Wi, but 
do not show the curves since they are very similar to those for Wq. 



6 Conclusion 

We have shown that the flux- form $ is the unique (ri— l)-form on a codimension-one saddle 
connection that describes the lowest order splitting of the manifolds upon perturbation. 
The integral of the one-way flux over a fundamental domain characterizes the transport 
across the manifolds in the perturbed system. For our example, the magnitude of the 
flux is strongly correlated with bifurcations in the homology of the crossing curves — near 
a bifurcation the flux is small. It would be nice to understand if this is a general feature 
of transport for volume-preserving maps. 
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Figure 10: Contours of the geometric flux through Wo as a function of ly and r for the perturbation (p8|). 
Also shown are bifurcation curves corresponding to the change in homology types of the zero contours of 
Mdj- In this case there is only one curve of bifurcation, corresponding to (1, 0) < — > (3, 1). 



In the future we also hope to study the evolution of the full manifolds numerically, 
to compare with our Melnikov results. We would also like to develop a nonperturbative 
method to compute the geometric flux, analogous to the action techniques for symplectic 
maps [1^. With this, we would like to verify that the geometric flux quantifies the 



transport observed numerically. 
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